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At the LHC Multiple Parton Interactions will represent an important feature of the mini- 
mum bias and of the underlying event and will give important contributions in many channels 
of interest for the search of new physics. Different numbers of multiple collision may con- 
tribute to the production of a given final state and one should expect important interference 
effects in the regime where different contributions have similar rates. We show, on the 
contrary, that, once multiple parton interactions are identified by their different topologies, 
terms with different numbers of multiple parton interactions do not interfere in the final 
cross section. 

PACS numbers: ll.80.La; 12.38.Bx; 13.85.Hd; 13.87.-a 

Keywords: Multiple scattering, Perturbative calculations, Inelastic scattering, Multiple production 
of jets 



1. INTRODUCTION 

The growing flux of partons at high energy will increase considerably the chances to have in- 
elastic events where more than a single pair of partons interact with large momentum exchange 
at the LHC[1J[2J|3J|4J|5J|6J|7J|8J. The phenomenon is originated by the increasingly large flux of 
partons at small fractional momenta. Once the final state is fixed the parton flux is maximized 
in the channel where the hard component of the interaction is maximally disconnected^]. In the 
resulting picture of the process, the dominant contribution at high energy is hence given by a 
set of Multiple Parton Interactions (MPI) where different pairs of partons collide independently 
in different points inside the overlap volume of the two interacting hadrons[10j[llj|12j. On the 
other hand, although the contribution with the largest number of initial state partons dominates 
at large hadron-hadron center of mass energes, a given final state may be generated by various 
competing processes, characterized by different numbers of partonic collisions [7]. Interactions 
with different numbers of partonic collisions populate in fact the final state phase space in a 
different way and one may always find kinematical regions where terms with different numbers 
of collisions give similar contributions to the cross section. In those kinematical regions, im- 
portant interference effects between different production mechanisms should be expected. The 
problem of interferences between terms with different numbers of collisions was, on the other 
hand, never discussed in the literature, while all theoretical estimates have always assumed inco- 
herence between contributions with different numbers of parton collisions, obtaining results not 
in contradiction with the available experimental evidence[13j[14J|15j|16j. 

The purpose of the present note is to gain some understanding of the problem by looking at 
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the kinematics of the different terms. After reminding the reader of the kinematical argument, 
which leads to the geometrical picture of MPI processes, we will analyze an interference diagram. 
The comparison between diagonal and off-diagonal terms in the cross section will allow to draw 
some general conclusions. 

2. MPI DIAGONAL SCATTERING DIAGRAM 

In Fig.l we show the cut diagram representing the contribution to the forward amplitude of 
a process with n partonic collisions. 




Figure 1: Unitarity diagram for the multi partem scattering cross section 



To study the kinematics we simplify the problem by limiting our discussion to the scalar case[T7]. 
We consider moreover all partons as identical particles. The soft vertices are assumed to be 
characterized by a non-perturbative scale of the order of the hadron radius R, which represents 
the (energy independent) scale of the transverse momenta and virtualities of the attached lines. 
To be definite, we consider the specific multi-parton interaction process where each elementary 
interaction T«, represented by the squares in the figure, generates two large pt partons with 
momenta pi and pi. 

Momentum conservation in the vertices limits the number of four dimensional variables to 3n — 1. 
Defining 

S% = ^(oi-Oi+i), Si = ^(ai-af i+1 ), Pi=pi+pi, (1) 
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in such a way that at + bi = Pi = a[ + b\, one may choose as independent variables: 



Pi, n — variables, 

Si, (n — 1) — variables, 

6'i, (n — 1) — variables, 

k or /, 



(2) 



since the overall momentum conservation is ^ Pi + k + 1 = Pa + -Pb- Sometimes the auxiliary 
dependent variables 



hbi-bi +1 ) = hPi-p i+1 )-5i 



(3) 



will be used. 

The four momentum variables will be represented in the light-cone form and the longitudinal 
and transverse components will be studied separately. The vertices <fi are "soft" and represent 
the non-perturbative partonic content of the hadron. In the hadron cm. one has 



1 

R 
1 

R 
1 



(4) 



a 



< 



where R is the hadron radius. By performing a boost to reach the cm. frame of the two 
interacting hadrons, one obtains 



a±, b ± < ± 

n h < 1 

"+> U ~ ~ 2M R 

_ u < 2M 1 



(5) 



where 5 = (Pa + -Pb) 2 is the hadron-hadron cm. energy and M is a scale of the order of the 
hadron mass. When looking at <5j_ and 5i+ one hence finds that both variables become smaller 
and smaller at large cm. energies: 



5- < 2Mi 

X. < 2M 1 



(6) 



The variables <5j_ are thus relevant only for the vertex 4>a and for the propagators of the lines 
with momenta a«, since the kinematical range of aj_ is of 0(2M/RVS). At the leading order, one 
needs in fact to take into account only of the kinematical variables which grow as VS in the hard 
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interaction vertices, while the lower vertex (pB and the propagators of the lines with momenta bi 
will remain practically constant for variations of <5_ of 0(2M/R\/S), as the kinematical range 
of bi- is of O(VS). Conversely the variables <5j+ are relevant only for the vertex <f>s and for the 
propagators of the lines with momenta bi, where all " + " components have a kinematical range 
of 0(2M/R\/S). Similar considerations hold for the variables a\, b\, b\_, S' i+ and for the vertices 
4>' A , <f>' B . One may hence define: 

n— 1 

, / i a f (j)A(ai . . . a n , k; j) yj dSi- 
f A {a 1+ ...a n+ ;a 1± ...a n± ;k-;j) = j - 2 11 ~2^T j. +=0 

, , , , -a f <t>B(bi ■ ■ -bnJ-J') j-r dS i+ 

il>B{bi-...b n -;b 1± ...b n± ;l+;j) = J j^-p [[ — _ =q (7) 

where k + , k± are given in terms of aj+ , ajj_ and /_ , l± are given in terms of bi- , bi± , while 
the value of fc_ is determined by the values of k + , k± and by the value of the invariant mass 
of the remnants of the hadron A, (Pa — k) 2 . l + is similarly determined by l± and by the 
value of the invariant mass of the remnants of the hadron B, (Pb — I) 2 - All other variables which 
characterize the remnants of A and B are labeled by the indices j and j' respectively. One has: 



ai + bi = Pi = a ■ + b'i 

a i+ + b i+ = a' i+ + b' i+ ~ a i+ ~ a' i+ = P i+ (8) 
cii- + bi- = a'i_ + b\_ ~ 6j_ ~ 6 ■_ = Pi- . 

The " + " and "-" components ^= are the "+" or "-" fractional momenta x^ 4 , xf and arc 
given by the final state observable quantities Pi + and Pj_. All longitudinal variables are hence 
either integrated or determined by the final state observables. In particular, taking into account 
only the terms which grow with y/S, one has that xf = x'/ 4 , xf = x'^ . 

As we are particularly interested in the structure of the interaction in transverse plane, we 
express the elementary interaction amplitude as a two-dimensional Fourier transform with respect 
to the relative transverse distance rf. 

T(\i,t i± ) = i- J e ^ • ^fiX^r^n (9) 

where U = \(ai — Pi — bi + pi) is the momentum transfer and all longitudinal variables are 
summarized by Aj. 

The multiparton cross section may hence be obtained by performing the n two-dimensional 
integrations on r^, r[ and in the following variables: 

Pi± = ail. + h± = a'i± + b' iL , n variables, 
* Qi± = ( a i± — bi±)/2, n variables, (10) 
q' i± = (a' i± - b' i± )/2, n variables. 
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Momentum conservation imposes however a constraint on the integration variables. As the 
incoming hadrons have vanishing transverse momenta one has 

k± + Y, a i± = k ± + Z~2 a '^ = > l ± + Y, bi± = Z -L + S 6 a = °> ( n ) 

i i i i 

which imply 

1>*-L " = ( 12 ) 

i 

in such a way that the cross section is given by 5n — 1 independent two-dimensional integrations. 
The integrals on the transverse components are suitably performed by representing the ip func- 
tions as two-dimensional Fourier transforms with respect to the transverse parton coordinates 

&i j S{, Sj , s j . 



ip A {xf ■ ■ ■ a 1± ... a n± ;(P A - k) 2 ;j) 

/i t 
i; A (xf . . . xt Sl . . . s n ; (P A - kf- j) [] e ia * 

^W-^; a'i± ■■■ a' n± ;{PA-k) 2 ;j) 

/i b 
r A (xt--xt s' 1 ...s' n ,(P A -k) 2 ;j)Y[e 

1> B (x?...x%; b 1± ... b n± ;(P B -l) 2 ;f) 

&,(sf ...x*;- Sl ... s n ; (P B - lf;f) J[ e ^ ' ^ 
P B (x?...x%; b' 1± ... b' n± ;(P B -l) 2 ;j') 



d 2 s- 

ldi± ■ Si a 6 » 

2tt 



■i i d 2 n' 
2vr 



i=l 



f...x^s[... s' n , (Ps - I) 2 ; f) e ~ ib i± " *i °£ 



2vr 

i=l 



(13) 



where the longitudinal components are expressed through the fractional momenta xf, x~ . The 
constraint in Eq.(12) is imposed by introducing the further integration 



_L_ J e iPJ:(Qi±-Qii_)d 2 (3 (14) 

where (3 is the hadronic impact parameter. One may integrate on n^J-> II dq'± by 

using the relations 

4± = \Pi± + Qi±, K± = \Pi± - «L ( 15 ) 

U± = qi± - \(Pi± -Pi±), t' i± = q' i± - ^(pij_ -pi±). 
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One obtains 



/ ^2 - S((s i + - Si - 4-^/2) 

S ^ ^si-si + n + p) 

5{s t i -s t i + T t i + (i) (16) 



(2vr) 2 

which is equivalent to <5(sj — — j(r- — r^)) — s'j + si — s' i )5(si — s~i + ri + (3). The integrations 
on r,; and r' involve the factor 



^ J e ^ ■ (*i-L " *i±)f (Ai^Of *(A i y i )rf 2 r i rf 2 rJ. (17) 

The difference \{r[ — ri) hence represents a measure of the localization of the interaction. As it 
appears from eq.(9), typical values of ri and r[ are of the order of 1/U±, which for large transverse 
momenta is equal to l/pi±. When the momentum exchanged in the elementary interaction is 
large, r, becomes much smaller as compared with the hadron size and it may be neglected, as 
well as rl, everywhere except in the expression in Eq.17. The dominant contribution at large cm. 
energy and at large exchanged momenta is thus obtained by making the following replacements 
in the evaluation of the discontinuity of the diagram in Fig.l: 



5(si - s'i - -(r- - ri))S(si - s- + s; - s-)<K -Si + n + 0) 

=> ${si- s'i) 5(si-$)6(si - Si + (3) 

j^yt J e^'^-^nX^i Xi ^d 2 nd 2 ^ 

=> iTix^x'ilpuPi)?, (18) 

where only the kinematical components which grow with 5 are taken into account in the eval- 
uation of T. One is hence left with the integrations on f cI(Pa — fe)_d(Ps — l)+d(3 f\ d,Si, while 
the inclusive cross section depends on |V>a| 2 an d \4>b\ 2 - Explicitly: 



an= 2^). j ^2,\i>A{xi ...x^;s 1 ...s n ;{P A -k) 2 ;j)\ 2 

3 

x J2\M X ? ■■■x*;s 1 -(3...s n -f3;(P B -l) 2 ;j')\ 2 
,/„ ,n ,/„ ,n d 2 8 r n 1 4 t-t d 2 Si 

x d(p A -ky d (p B -i) +J ^- 2 [^} 

x \T(xf,x?\ P i,Pi)\ 2 d$i dx?...dx£dx? ...dxl (19) 
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where $j is the invariant adimensional final state phase space of the elementary interaction 
CLi + hi — > pi +pi and all elementary interactions are considered as indistinguishable. The factors 
2, 7r etc. originate from the Jacobian which leads to the variables in Eq.19. 

After multiplying and dividing by the flux factors 2Sxfxf , one may introduce the elementary 
partonic cross sections 

a(xf,xf;p cut ) = — / \T(x?,x?\p t ,p t )\ 2 d<S> t (20) 

where p cu t is a cutoff in the transverse momenta of final state partons, introduced to allow to 
compute the cross section in perturbative QCD. The multi-parton densities T(xf, Sj) are hence 
defined as 



r , .In 2 Sn-iflZxi 

T(x 1 ...x n ;s 1 ...s n ) - 



(27r) n+l 4 „-l y/2(l-^Xi) 



I Y J \^--^n,S 1 ...S n ;{P-k) 2 ;j)\ 2 d{P-k) 2 (21) 
J 3 

and the cross section is finally expressed by 



a n = — T A (x 1 ...x n ;s 1 ...s n ) 

n\ J 

x l[a(xf,xf;p cu t)T B (xf ...x%;si - [3 . . . s n - [3)d 2 [3d 2 Si dxfdxf , (22) 

i 

which represents the superposition of n elementary collisions localized in regions with transverse 
size of the order of l/p cu t, much smaller as compared with the hadron size, and with the mean 
value of the transverse coordinates Sj. 



3. INTERFERENCE TERMS 



There are two different ways of producing interference terms. A possibility is to have a 
different number of hard collisions on the left and on the right hand side of the cut. Another 
possibility is to have the same number of hard collisions on both sides of the cut, in which case 
interferences are produced by reshuffling the final states of the hard collision terms. 

We will first consider a case of interference between a term with n and a term with n — 1 
collisions, the corresponding unitarity diagram is shown in Fig.2. 
In this case the independent variables are 3n — 2: 



Pi n — variables, 

Si (n — 2) — variables, , N 

V ' ' (23) 

5- (n — 1) — variables, 

k or I. 
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The longitudinal variables may be discussed as in the previous case. The obvious difference is in 
the resulting expression, which is no longer a modulus squared. The main point for the present 
considerations concerns the integrations on the transverse variables. Analogously to the diagonal 
case, the integrations on the transverse variables are conveniently discussed taking the Fourier 
transforms of the functions tjj: 



?pA,n-i {x , xi ... x n ;a _L,a 3 _L . . .a nJ _; (Pa - k) ;j) 



/n ,2 
4>A,n-i(x$,i4 ...afi; so,s 3 ...s n ;(P A -k) 2 ;j)]J e iai± ' Si - - 



i=l 

ip B ,n-i{xo,Xz ... x^; b ±,b 3± . . .b n± ; (P B - Z) 2 ;j) 

= [$ Bjn „ 1 (x B ,x$...xZ ] -s ,S3...s n ;(PB-l) 2 ;j)fle ib ^-^,(24) 

where a label representing the number of interacting parton lines (n — 1 in this case) has been 
introduced. The complex conjugate functions ^* are the same as in the diagonal case. As 
discussed in the diagonal case, the large transverse momenta exchanged in each elementary 
collision localize the interactions in transverse space regions much smaller as compared to the 
hadron radius. While the discussion of the transverse variables may be done following the lines 
of the diagonal case, the treatment may be simplified by neglecting since start the distances 
between the interacting partons n, r\ in comparison with the parton coordinates Si, Sj, s^. 
Introducing the variable qo = ao — bo one has 



2tt 
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'a =l(P 1 +P 2 + q ) 
b = \(P l +P 2 -q ). 



The integration on qo± gives 



J dq 0± -► S(s - s + 0), (26) 
the integrations on q[, q' 2 (previously defined) give 

J dq ' 1± ^ s^-^+P) 

dq' 2 ± 5(s' 2 - 4 + P) (27) 



and all other integrations on <^ are the same as in the diagonal case. At this stage the functions 
ip depend on the transverse variables as follows (the dependence on the fractional momenta and 
on the invariant mass of the remnants of the hadron is implicit) 

Vu,n-l(so,S 3 • ■■M,n( S 'l' S 2) s 3 ■ • • )^B,n-l( s ~ P,S 3 ~ P ■ • • )4>*B,n( s l ~ P, 4 ~ P, 4 ~ P ■ ■ ■ )■ 

(28) 

One may now perform the integrations on P±±, P 2 ±. The result is 



J dPi± -> <S(s + so - 4 - 4) 
y <^P 2 _l -»■ <5(s + s - 4 - 4)> 



(29) 



which is equivalent to 5(sq — s' 1 )5(sq — s' 2 ). All other integrations on give the same result as 
in the diagonal case. One hence obtains that two transverse variables coincide in ij)* A n and ip* B n 
and the cross section is proportional to the integral 



/ 



tjj A ,n-l(s , S 3 ... S n )i)\ n (so, S , S 3 . . . S n ) 

n 

1pB,n-l{so - P,S 3 - P...S n - P)i>* B , n( s ~ P, SQ ~ P, S 3 - P . . . S n - P)ds dP JJ ds 



i 

(30) 



Analogously to the diagonal case, sq, s 3 . . . s n represent the transverse coordinates of the positions 
of the interaction regions. One may hence conclude that, in the interference term between a n 
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and a n — 1 collisions amplitude, the hard component of the interaction is localized in n — 1 
points in transverse space. 

When the number of hard collisions is the same in both sides, a non-diagonal contribution 
may be obtained by linking, across the cut in Fig.l, two different collision amplitudes through 
the produced large pt final states. In such a case one obtains that the positions of the two hard 
interactions are localized within the same region, with size of 0(r±), which implies that the 
number of integrations in the transverse coordinates is reduced by one unit with respect to 
the diagonal case. Also in this case the interference term hence corresponds to a case where 
the hard interaction is localized in n — 1 points in transverse space. Rather obviously further 
crossings of the hard (p, p) lines would further reduce the number of transverse integrations and 
hence the number of points in transverse space where the hard component of the interaction is 
localized. 

4. CONCLUDING DISCUSSION 

A given multi-partons final state may be produced by interactions involving different numbers 
of partons in the initial state and the cross section, resulting from the coherent sum of all different 
terms, is expressed by a sum of diagonal and off diagonal contributions. As shown in Sec. 2, the 
diagonal contribution, corresponding to a term with n partons in the initial state, is given by the 
incoherent super-positions of n disconnected parton interactions, localized in n different points 
in transverse space. As shown in Sec. 3, the hard component of the interaction corresponding to 
off diagonal contributions is disconnected and localized in no more than n—1 points in transverse 
space. 

One may hence argue that interference terms do not represent corrections to the n-partons 
scattering inclusive cross section. They rather correct the n — 1-partons (or less) scattering 
inclusive cross section. 

Partons are in fact localized in the hadron by the momenta exchanged in the interaction. 
When partons are localized inside non overlapping regions, much smaller as compared to the 
hadron size, they are only connected one with another through soft exchanges and the picture 
of independent parallel collisions described in section 2 is a meaningful one. If, on the contrary, 
partons are localized by the interaction inside overlapping regions, much smaller as compared to 
the hadron size, they are allowed to interact by exchanging momenta of the size of their virtuality, 
which implies that the evaluation of the interference term, as discussed in section 3, is no longer 
adequate. 

In Fig. 3 we show an interference diagram between a single collision amplitude, where two 
partons interact at tree level producing four large pt partons, and an amplitude, where two 
partons, generated by the same short distance quantum fluctuation in the hadron B, interact 
with two partons of the hadron A, producing two pairs of large pt partons. 
As it appears looking at Fig. 3, because of the localization of the hard component of the interac- 
tion, the problem of interference is strictly linked to the problem to evaluate the single scattering 
amplitude, at higher orders in the coupling constant and including higher twists in the hadron 
structure. 

One may hence conclude that a convenient way to distinguish the different terms in a MPI 



11 




Figure 3: A particular case of interference diagram 



process is by their different topologies: 

As a consequence of the different scales involved in the interaction, namely the hadron size and 
the large momenta exchanged, the structure of the hard component may be disconnected, with 
the different hard parts linked only through soft exchanges. The disconnected parts of the hard 
interaction are localized in different regions in transverse space, inside the overlap of the matter 
distribution of the two interacting hadrons. The different MPI terms are to be understood as the 
contributions to the final state due to the different disconnected parts of the hard component of 
the interaction. 

In the instance of a single scattering the hard component is wholly connected and hence 
localized inside a single region in transverse space. In the simplest case, a single interaction 
is well described by the simple QCD-parton model recipe, with the parton interaction cross 
section evaluated at the lowest order in the coupling constant and convoluted with the parton 
distributions. In other cases, one may need to take into account higher order terms in the 
coupling constant to evaluate the partonic cross section or/and to include higher twist terms in 
the parton distributions. The effects of higher order terms and of higher twists may be more 
important when dealing with final states with several large pt partons. 

A multiple interaction, on the other hand, has to be understood as a process where the hard 
part of the interaction is disconnected and localized in a number of different regions in transverse 
space. In the simplest case, in each different region the interaction may be evaluated at the lowest 
order in the coupling constant. 

The main observation in the present note is that when MPI are understood in the topolog- 
ical sense described here above, different MPI terms, corresponding to different localizations in 
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transverse space, namely to different topologies, do not interfere and the final cross section is ob- 
tained simply by the superposition of the cross sections, due to the contributions of the different 
topologies of the hard component of the interaction. 

The topological feature represents on the other hand also the property which allows to recog- 
nize the contribution to the final state due to MPI. In each single parton collision all transverse 
momenta need to balance and a MPI process contributes to the cross section generating different 
groups of final state partons where the large transverse momenta are compensated separately. 
The compensation of transverse momenta within different subsets of large pt final state partons 
was in fact the feature which allowed the experimental identification and study of double parton 
collisions[13-16]. Interestingly, in a recent study of MPI at LHC energies it was shown that, 
requiring the compensation of transverse momenta within different subsets of large pt final state 
partons, one my expect to be able to isolate contributions due to triple parton collisions also in 
channels with a relatively low cross sections and with relatively large transverse momenta[7]. 
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